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ABSTRACT 

B2 0954+25A, detected by the Fermi satellite, is a blazar with interesting observational 
properties: it has been observed to transit from a jet dominated to a disk dominated 
state; its radio spectrum appears flat at all observing frequencies (down to 74 MHz); 
optically, the H/3 line profile is asymmetric. The flatness of radio spectrum suggests 
that the isotropic emission from radio lobes is very weak, despite the large size of its jet 
(> 500 kpc). Its broad-band spectral energy distribution is surprisingly similar to that 
of the prototypical 7-ray, radio loud. Narrow Line Seyfert 1 (7-NLSI) galaxy PMN 
J0948+0022. In this work we revisit the mass estimates of B2 0954+25A considering 
only the symmetric component of the H/3 line and find (1-3) xlO^ Mq. In light of our 
composite analysis, we propose to classify the source as a transition object between 
the class of Flat Spectrum Radio Quasar and 7-ray, radio loud NLSl. A comparison 
with two members of each class (3C 273 and PMN J0948+0022) is discussed. 

Key words: quasars: individual: B2 0954+25A- galaxies: jets - 7-rays: observations. 



1 INTRODUCTION 

A major breakthrough in the comprehension of the AGN 
phenomenon lias been the formalization of the so-called 
"unified model" (Urry & Padovani 1995). AH active galactic 
nuclei are now believed to be similar objects, whose gross 
observational features depend on the black hole mass, the 
accretion rate, the viewing angle and the presence of a jet. 
Our understanding of finer details is, however, far from com- 
plete: each source has its own list of peculiarities that are 
difficult to explain coherently in the framework of the uni- 
fied model. Furthermore, the values of the black hole mass, 
that can be inferred only indirectly, are still uncertain. Iden- 
tification and analysis of peculiar features in single sources 
is thus important to improve our understanding of the AGN 
phenomenon. In this context the Fermi blazar B2 0954-I-25A 
(= OK 290 = TXS 0953-^254), detected in the 7-ray band 
by the Fermi satellite, has several interesting features worth 
investigating. Interest for this source has been stimulated by 
the wide range of values for its black hole mass and width 
of hydrogen Balmer lines found in literature. 

At optical wavelengths the source has been observed in 
at least three emission states: the partially jet dominated 
state in 1987 (Jackson & Browne 1991), as observed by the 
Isaac Newton Telescope (INT); the disk dominated state in 
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2004, as observed by the Sloan Digital Sky Survey (SDSS, 
Abazajian et al. 2009); and the completely jet dominated 
state in 2006, observed again with the SDSS. The optical 
spectrum reveals that the H/3 emission line shows an asym- 
metry in its red wing, as observed in many AGNs. This 
asymmetry is one of the likely reasons for the different val- 
ues of the black hole mass found in the literature. Depending 
on the decomposition made to derive the real width of the 
line, this object could be classified as a Narrow Line Seyfert 
1 (NLSl) characterized by a full width at half maximum 
(FWHM) smaller than 2000 km s'^ (Goodrich 1989). In 
this case the source would join the small group of radio-loud 
NLSl that have been detected by the Fermi satellite (Abdo 
et al. 2009a,c; Foschini 2011). Indeed, its spectral energy 
distribution (SED) is almost identical to PMN J09484-0022 
(Abdo et al. 2009b; Foschini et al. 2009), the first radio loud 
NLSl detected by Fermi. 

In this work we present a multi-wavelength study on B2 
0954-I-25A, with particular emphasis on the spectral analy- 
sis at optical wavelengths, broad-band SED modeling and 
radio properties. Our aim is to provide a coherent picture 
of its physical properties, according to observational data. 
Our data analysis relies on data from several facilities, as 
discussed in the following sections. General data are shown 
in Tab. 1. 

Throughout the paper, we assume a ACDM cosmology 
with Ho = 71 km s"^ Mpc"\ f^m = 0.27, = 0.73. 
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Table 1. General data. 



Parameter 


Value 


Units 


R.A. 


149.208 


deg 


Declination 


25.254 


deg 




0.70747 




dL 


4.3 


Gpc 


Ab^ 


0.158 






?, 


102" 


E(B \-)<' 


0.0:-!77 




From [OIII] line fitting (sec 





Table 2. Observation dates for all data used to build the SED in 
Fig. 2. 



^ Galactic absorption from Schlegel et al. (1998) 

Computed using the LAB survey (Kalberla et al. 2005) 
From http : //galex . stsci . edu 



2 THE SOURCE B2 0954+25A 

The source B2 0954+25A (a=0.712, Burbidge & Strittmat- 
ter 1972) is a compact, radio-loud, flat-spectrum radio 
quasax (FSRQ). It has been frequently used in statistical 
studies on radio properties of quasars since it is a relatively 
luminous radio source (>1 Jy, Kuchr ct al. 1981), whose 
emission extends to very low radio-frequencies (74 MHz, 
Cohen et al. 2007). The radio spectrum is usually flat and 
becomes inverted during burst activity (Torniainen et al. 
2005). The jet is clearly visible in several radio maps: see e.g. 
the VLA radio maps at 1.64 GHz in Murphy et al. (1993) and 
the VLBA radio maps at 22 and 43 GHz in Lister & Smith 
(2000). The core component has angular size ~ 0.23 x 0.07 
mas (at 15 GHz), corresponding to a linear size of 1.6 x 0.49 
pc (Kovalev et al. 2005) . Several components in the jet show 
superluminal motion (up to 12c, Kellermann et al. 2004). A 
one-sided jet (projected size ~50 kpc) extends from the core 
in the south-west direction (Liu & Zhang 2002). 

In the optical band, the source is unresolved, vari- 
able (Pica et al. 1988) and slightly polarized (1.29%, Wills 
et al. 1992), with m, ~18 mag. The bolometric luminos- 
ity (estimated from SED fitting, Woo & Urry 2002) is 
log(I/boi/erg s~^) = 46.59. Virial black hole mass estimates 
arc log(M/MQ)=8.7 (Liu et al. 2006) and log(Af/M0)=9.5 
(Gu et al. 2001). Both these estimates rely on a FWHM es- 
timate of the H|0 emission line given in Jackson & Browne 
(1991), who found FWHM(H/3) = 65 A (rest frame), cor- 
responding to ~4000km s^^. The source is also present in 
the Shen et al. (2010, hereafter SIO) catalog, which contains 
several measures obtained by automatically analyzing the 
SDSS/DR7 spectrum. SIO report a FWHM(H/3) = 1870 km 
s"^ and log(M/MQ)=8.6 (computed with H/3), or 9.3 (com- 
puted with Mgll). The spectrum analyzed in SIO is exactly 
the same as the one we use in §4. 

Several X-ray facilities {Swift/XKT, Chandra, ROSAT, 
Einstein) observed the source at different times measuring 
fluxes in the range (2.5-20) xlO"^^ erg cm~^ s~^ (§3 and 
Fig. 2). Finally, B2 0954-I-25A is present in both the lyr 
and 2yr Fermi Large Area Telescope (LAT) point source 
catalogs (Abdo et al. 2010, 2011, with catalog names IFGL 
J0956.9-I-2513 and 2FGL J0956.9-I-2516 respectively). 



Facility 



Einstein 

INT 

2MASS 

SDSS (photometry) 

SDSS (spectrum) 

GALEX 

Swift 

Chandra 

Fermi (average) 



Date of obs, 
1979- 
1987 
1998 
2004 
2006- 
2006 
2007 
2009 
Aug 2008 



11- 06 

12- 16 

11- 30 

12- 13 

01- 05 

02- 05 
06-01 
01-20 
- Jun 2011 



3 BROAD-BAND SED 

To build the broad-band SED (Fig. 2) we collected pub- 
licly available data from several facilities using NED.^ Grey 
crosses are from radio facilities, black squares are from Sloan 
photometry, red butterfly is from ROSAT, purple diamond 
is from Einstein. We added the data from i^ermi/Large Area 
Telescope (LAT) (green and black butterflies), Swift/ji- 
Ray Telescope (XRT) (blue dots), Chandra (green butter- 
fly), GALEX (orange open circles), 5uii/t/Ultra Violet Opti- 
cal Telescope (UVOT) (blue circles) and 2MASS (red tri- 
angles). Tab. 2 shows the observation dates for all facilities. 
Here we provide a brief description of data analysis for each 
of the facilities: 



Fermi/LAT — The 2LAC catalog (Ackermann et al. 2011) 
reports a flux for the source of (2.1 ± 0.2) x lO"'^^ erg 
cm"^ at 618 MeV, and a photon index F = -2.39 ± 0.08. 
We also ran a likelihood analysis on the first 34 months 
of data collected by Fermt/LAT (using ScienceTools ver. 
9.23.1 and Instrument Response Function P6_V3) with time 
bins of one month, to check for the occurrence of major 
flares. We found no significant variability of the flux nor of 
the spectral index, i.e. a constant flux and spectral index 
model fits well observational data and cannot be rejected on 
statistical basis. The flux in the 0.1-100 GeV band averaged 
over 34 months is compatible with the value reported in 
2LAC. 

Swift/XRT — On June 1st, 2007 Swift performed a pointed 
observation (obsID 00036325002) of a nearby source, thus it 
has been possible to extract data for B2 0954-I-25A. We ex- 
tracted a spectrum using the xrtpipeline script and binned 
it to have at least 20 counts in each bin. Then we used 
XSPEC to fit a simple power law with the hydrogen column 
fixed to the Galactic value Nn = 3.2 x 10^° cm"^. The re- 
sulting de-absorbed flux in the 0.3-10 keV energy band is 



(2.5 ± 0.3) X 10"^^ 
r = 1.74 ± 0.09. 



erg cm ^ s ^and the photon index is 



Chandra — On June 20th, 2009 Chandra pointed our source. 
We extracted a spectrum using the specextract script, with 
a minimum of 30 counts in each bin, then performed a fit 
against a simple power law using Sherpa. Again the A^h 
column has been kept fixed to 3.2 x 10^" cm"'^. The re- 
sulting de-absorbed fiux in the 0.3-10 keV energy band is 



^ Nasa/IPAC Extragalactic Database http : //http : / /ned . ipac . 
caltech.edu/ 
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(1.32 ± 0.06) X 10 erg cm ^ s ^and the photon index is 
r = 1.75 ± 0.07. 

GALEX — on Feb 5th, 2006 GALEX observed the source 
and provided the following photometric measurements: 
(269 ± 18) /iJy at 1528 A (far UV) and (535 ± 14) ^ijy 
at 2271 A (near UV). Using a colour excess E(B-V) = 

0. 0377 we computed the de-reddened fluxes using CCM 1989 
parametrization (Cardelli et al. 1989). The resulting uF„ 
values are (7.00 ± 0.47) x 10"^^ erg cm"^ s"^(FUV) and 
(9.71 ± 0.26) X 10-^2 erg cm^^ s-^(NUV). 

Swift/UVOT — - An aperture photometry analysis of the u 
filter data provided a count rate of 3.16 ±0.11 s~^. The cor- 
responding (Poole et al. 2008) de-reddened (Cardelli et al. 
1989) vF^ flux at 3501 A is (2.14 ± 0.07) x 10"^^ erg cm"^ 

2MASS — An IR observation of the source has been per- 
formed on Nov 30th, 1998 and results are available in 
the 2MASS catalog. Magnitudes are mj = 16.53 ± 0.13, 
■mn = 16. 12 ±0.19 and ttik = 15.34 ±0.14. The correspond- 
ing ^ uFv fluxes are: (9.5 ± 1.1) x 10~^^ erg cm~^ s~^at a 
1.2 ^im (J band), (6.6 ± 1.1) x 10"" erg cm"^ s"^at a 1.6 
^im (H band) and (6.8 ± 0.8) x 10"" erg cm"^ s"^at a 2.2 
/im (K band). 

Finally, we also added the spectroscopic measurement 
from Sloan/SDSS (orange lines, see also §4) taken about one 
year after the SDSS photometry and about one month before 
GALEX observations. Looking at the SED in optical/UV 
it is clear that the source raised its flux by a factor ~10 
in about one year. In June 2007 the source returned to a 
lower state (in optical/UV) as shown by the Swift/WOT 
measurements. Thus, the source shows at least a "low" and 
a "high" state in optical/UV. A similar trend, although of 
smaller magnitude, is shown at X-ray wavelengths with the 
lowest flux being measured by ROSAT and the highest being 
measured by Einstein. 

3.1 Radio properties 

B2 0954+25A shows a flat radio spectrum from ~ 10* Hz to 
~ 10^^ Hz. At 5 GHz, the luminosity of the core is log(Z/c/erg 
s"^ Hz"^) = 34.5, while that of the extended region is 
log(Le/erg s"^ Hz"^) = 32.4, thus the source is highly core- 
dominated with a core dominance parameter 7?cd ~ 100, 

1. e. at the higher end of the core-dominance parameter dis- 
tribution (Browne & Murphy 1987). The weakness of the 
extended emission is also evident at low radio frequencies, 
because the radio spectrum remains flat down to 74 MHz, 
with a spectral index a ~ —0.2 between 74 and 365 MHz 
{Fi, Qc v°'). Assuming a spectral index for the optically thin 
emission in the extended region of -0.7, and extrapolating 
the flux from 5 GHz, the frequency at which the isotropic 
and the beamed jet luminosities are equal is ~10 MHz, in 
agreement with the observed overall flat radio spectral in- 
dex. This suggests that the jet is well aligned along the line 
of sight. 

^ http : / /www. ipac . caltech.edu/2mass/releases/allsky/faq. 
html#j ansky 
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Sp, Index [76 MHz, 365 MHz] 

Figure 1. Distribution of spectral index between the (observed) 
frequencies 74 and 365 MHz for three different blazar catalogs: 
CRATES (Healey et al. 2007, blue dashed line), BZCAT (Mas- 
saro et al. 2009, red dotted line) and Fermi 2yr catalog of AGN 
(Ackermann et al. 2011, 2LAC, black solid line). Thicker lines 
show cumulative fraction (values on right axis). Radio fluxes to 
compute the spectral index are taken from the VLA survey at 
74 MHz (Cohen et al. 2007) and the Texas survey at 365 MHz 
(Douglas et al. 1996). 



We can compare the 74-365 MHz radio spectral index a 
of B2 0954±25A with the distribution of the same spectral 
index of other blazars. To this aim we show in Fig. 1 the 
distribution of spectral index between the (observed) fre- 
quencies 74 and 365 MHz for three different blazar catalogs: 
CRATES (Healey et al. 2007, blue line), BZCAT (Massaro 
et al. 2009, red line) and Fermi 2yr source catalog of AGN 
(2LAC, ah "CLEAN" sources, Ackermann et al. 2011, black 
line). The radio fluxes to compute the spectral index are 
taken from the VLA survey at 74 MHz (Cohen et al. 2007) 
and the Texas survey at 365 MHz (Douglas et al. 1996). 
The number of sources detected at both frequencies is 1282 
for the CRATES catalog; 605 for the BZCAT and 218 for 
the Fermi 2LAC catalog. The distribution shows that only 
~30% of the highly beamed blazars detected by Fermi have 
a steep (q < —0.5) spectral index between 74 and 365 MHz. 
Radio-selected blazars from CRATES show a broader dis- 
tribution with ~50% sources having spectral index <-0.5. 
The BZCAT distribution is similar to CRATES, although 
with a larger fraction of very flat or inverted radio sources. 

The extended region of B2 0954±25A, as observed in 
the 1.64 GHz VLA radio map, has a projected size of ~50 
kpc (Murphy et al. 1993; Liu & Zhang 2002). This is rather 
large, especially considering that the viewing angle 9^ is 
small, as suggested by the observed superluminal motion 
(12c, Kellermann et al. 2004) and by the SED modehng 
(§3.5 and Tab. 3). Assuming 9^ in the range 3-6 degrees, 
the de-projected size is then in the range 0.5-1 Mpc, typi- 
cally the size of a giant radio lobe. 

We conclude that the radio properties of B2 0954+25A 
indicate a highly core dominated source, with an extended 
component that, albeit reaching the dimension of a giant 
radio lobe, is rather weak, making the core dominance pa- 
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Figure 2. Broad band spectral energy distribution, and related models, of B2 0954+25A. Data are provided by several facilities (see 
§3). Upper panel shows all available data, from radio to 7— ray energies, as well as the SED models in the low and high synchrotron state. 
The lower panel shows a detailed view of the SED from IR to 7-rays, and highlights the most important components of the models. 



rameter very large, at the high end of the distribution values 
for blazars. 



3.2 Modeling the SED 

To model the SED we use the one-zone leptonic model de- 
scribed in detail in Ghisellini & Tavecchio (2009), whose 
main characteristics are briefly described in the Appendix. 
The SED of the source is shown in Fig 2. The lower panel 
shows a zoom of the optical. X-ray and 7-ray data part of 
the SED, together with the models used to interpret the two 
states of the source. As discussed above, SDSS observed B2 
0954-I-25A in two occasions, approximately one year apart. 
In the first observation the SDSS made the photometry in 
its usual optical filters (black squares in Fig. 2), while the 



spectrum, one year later, revealed a very high state, dom- 
inated by the beamed synchrotron emission (orange lines), 
in agreement with the GALEX observations (orange open 
circles) . 

For the low optical state we have simultaneous UVOT 
and XRT Swift observations: the UVOT point (blue fiUed 
circle) agrees with the photometric SDSS data points (black 
squares) and with the 2MASS data (red triangles). We have 
then modeled these data together with the XRT X-ray data, 
assuming also that the average Fermi/LAT flux and spec- 
trum is a good representation of the high energy emission of 
this state. 

On the other hand we have no information about the 
level of the 7-ray flux for the high optical state of the source. 
We therefore assume that it corresponds also to a high state 
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of the X ray emission, that can be represented by the max- 
imum X-ray flux observed (violet diamond point), and to 
the same 7-ray flux. 

The Chandra and Fermi data shown as green ties are 
simultaneous (even if the Fermi spectrum is the average of 
one week of observations), but we do not have any informa/- 
tion about the optical-UV state. Therefore we did not try 
to model these data. 

In Tab. 3 we list the relevant parameters used to inter- 
pret the two states of the source, namely the disk and the jet 
dominated one. Tab. 4 reports the jet powers corresponding 
to the two states, calculated as explained in the Appendix. 



State 


logPr 


logPB 


logPe 


logPp 


jet dominated 


44.81 


45.75 


44.92 


47.71 


disk dominated 


44.56 


44.77 


44.71 


47.41 


PMN J0948+0022 


45.30 


44.35 


44.71 


46.68 


3C 273 


45.05 


46.09 


44.90 


47.48 



Table 4. Jet power in the form of radiation, Poynting flux, bulk 
motion of electrons and protons (assuming one cold proton per 
emitting electron) for of B2 0954+25A (first two lines), PMN 
J0948-I-0022 and 3C 273 (Ghisellini et al. 2010). 



3.2.1 Guidelines for the choice of the parameters 

Accretion disk — The luminosity of the accretion disk can 
be estimated through the luminosity of the broad emission 

lines Lblr, assuming a given value for the covering factor. 
As will be discussed in detail in the next section, we can 
estimate Lblr ~ 6.3 x 10** erg s~^. A.s.suming a covering 
factor of ~0.1, the disk luminosity is of the order of 10*'' erg 
s~^. Given the uncertainties about Lblr, derived using only 
the Hp line, the covering factor and the mass of the black 
hole (controlling the La /Lion ratio) , this estimate necessarily 
has an uncertainty of at least a factor 2. 

Location of the emitting region — The SED in the low state 

shows that the high energy hump dominates the bolometric 
output, as in other powerful Flat Spectrum Radio Quasars 
(FSRQs). In these objects it is likely that the emitting region 
is located inside the BLR, since in these cases the seed pho- 
tons produced by the BLR, whose flux is seen enhanced in 
the comoving frame of the emitting region, dominate the ra- 
diative cooling. We then assume i?diss < -Rblr, where i?diss 
is the distance of the emitting region from the black hole, 
and ilsLR ~ lO^^Ly^g cm is the size of the BLR. 

Magnetic field — The value of the magnetic field is chosen 
in order to reproduce the ratio between the synchrotron and 
the inverse Compton flux. 

Bulk Lorentz factor — This has been chosen in order to 
reproduce the superluminal apparent velocity (,0app ~ 12, 
Kellermann et al. 2004). This also fixes the viewing angle, 
that must be of the order of l/F. Note that the value of F 
determines the value of the radiation energy density of the 
external seed photons (oc F^) and hence the value of the 
magnetic field required to have the observed synchrotron to 
inverse Compton luminosity ratio. 

Injected power — The power is injected throughout the 
source in the form of relativistic electrons. Through the con- 
tinuity equation we calculate the particle distribution as a 
result of injection, cooling and possible pair production. The 
total injected power is such that the radiation produced by 
these particles agrees with the observed data. The injected 
distribution is assumed to be a smoothly broken power law. 
The resulting distribution, modified by cooling, must agree 
with the observed slopes. 



3.3 Results 

The most notable feature of B2 0954-f25A is the large vari- 
ability in the optical, from a disk to a jet dominated emis- 
sion. This is interpreted as a change in the magnetic field, 
from ~ 5 to 20 G, accompanied by an increase of the total 
injected power (by a factor 2), and a decrease in the max- 
imum particles energy and a flattening of the high energy 
slope of the particle distribution. The location of the emit- 
ting region is somewhat less in the high optical state. We 
maintained the same bulk Lorentz factor. The power of the 
jet of B2 0954-I-25A carried by the magnetic fleld is larger in 
the high optical state (factor 10), while the differences in the 
other forms of power are minor. If there is one (cold) proton 
per electron, the power is dominated by the bulk motion of 
protons, and is of the order of (2-5) x 10*^ erg . The min- 
imum value of jet power is the one needed to produce the 
observed radiation, that is of the order of Pr ~ LobB/r^ = 
(3-7) X 10** erg s'^ 

These values of the physical parameters are similar to 
those found by fitting Fermi blazars (Ghisellini et al. 2010), 
with the black hole mass being in the low end of the distribu- 
tion. We have used a black hole mass M = 1.5 x 10*Mq. We 
have also fitted the spectrum with larger black hole masses, 
and obtain equally good fits with M = 3 x lO^M© (and 
Ld/LEdd = 0.34) and M = 5 x 10* Mq (and Ld/Lfidd = 
0.17). Both these alternative solutions are acceptable since 
they provide a sufficient Lion to photoionize the BLR. If 
M = lO^M©, instead, the computed disk emission shifts to 
smaller frequencies (i.e. the disk becomes colder) and then 
we are obliged to decrease its overall luminosity. As a con- 
sequence, we cannot reproduce the Swift/UYOT u point, 
and the ionizing luminosity becomes too small to efficiently 
photoionize the BLR. 

The picture that emerges is that of a typical Fermi 
blazar, with a fiat radio spectrum produced by a powerful 
jet, and black hole mass significantly smaller than 10^ M©. 
The source periodically cycles through different emission 
states, characterized by different amount of power emitted 
by the jet. Although the SED parameters are only rough es- 
timates of physical quantities, the possibility to model two 
different states of the same source provides at least a clue 
that the parameters kept fixed in both models (black hole 
mass, disk luminosity, BLR radius, covering factor) are reli- 
able. 



© 2011 RAS, MNRAS 000, 1-14 



6 G. Calderone et al. 



State 


J KllSS 


M 






p.' 


B 


r 








-^1 


^2 


O'cool 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 




[12] 


[13] 


[14] 


jet dominated 


45 (1000) 


1.5e8 


11.3 (0.58) 


335 


0.014 


21 


13 


3 


100 


2.5e3 


2 


2.5 


7 


disk dominated 


59 (1300) 


1.5e8 


11.3 (0.58) 


335 


7e-3 


5.2 


13 


3 


100 


5e3 


2 


2.9 


21 


PMN J0948+0022 


2 (1600) 


1.5e8 


9 (0.4) 


300 


0.024 


3.4 


10 


6 


800 


1.6e3 


1 


2.2 




3C 273 


120 (500) 


8e8 


48 (0.4) 


693 


0.015 


11.6 


12.9 


3 


40 


2e4 


1 


3.4 





Table 3. Input parameters used to model the SED of B2 0954+25A (first two lines), PMN J0948+0022 and 3C 273 (Ghisellini et al. 
2010). Col. [1]: State of the source; Col. [2]: distance of the blob from the bla<;k hole in units of 10^® cm and in units of the Schwarzschild 
radius; Col. [3]: black hole mass in solar masses; Col. [4]: Disk luminosity in units of 10*® erg and in Eddington units; Col. [5]: 
distance of the BLR in units of 10^® cm; Col. [6]: power injected in the blob calculated in the comoving frame, in units of 10*® erg s~^; 
Col. [7]: magnetic field in Gauss; Col. [8]: bulk Lorciitz; Col. [9]: viewing angle in degrees; Col. [10], [11] and: minimum and maximum 
random Lorentz factors of the injected electrons; Col. [12] and [13]: slopes of the injected electron distribution [Q(7)] below and above 
7b; Col. [14] minimum random Lorentz factor of the electrons cooling in one light crossing time R/c. The disk has an X— ray corona of 
luminosity Lx = 0.3L(j. The spectral shape of the corona is assumed to be ex exp(— /iv/lSO keV). 



4 OPTICAL SPECTROSCOPY 

We used optical spectra from the Isaac Newton Telescope 
(INT) and from SDSS^ (see Fig. 3). The spectrum from INT 
(observed on 16th Dec. 1987) has been derived directly from 
the plot given in Jackson & Browne (1991, their Fig. 2), thus 
it is suitable only for a qualitative analysis. The spectrum 
from SDSS (observed on 5th Jan. 2006) has been automati- 
cally flux- and wavelength-calibrated by the SDSS pipeline 
(Stoughton et al. 2002), de-reddenod using a Milky Way- 
like extinction curve (Pei 1992) and using Ab = 0.158 mag 
(Schlegel et al. 1998), transformed to rest frame of reference 
using our redshift estimate (^=0. 70747, §4.1) and rebinned 
by a factor of 2 (Fig. 3, see Fig. 5 for a closer view of the 
H/3 and [OIII] region). 

Fig. 4 shows the profiles of the main emission lines, 
namely Mgll (2800 A), (4863A) and [OIII] (5007A), 
in velocity space. The spectrum has been continuum- 
subtracted and all profiles have been normalized to have 
maximum of 1. Finally all profiles have been smoothed by a 
~4 A boxcar average to provide a clearer view. The dashed 
lines allow to approximately measure the full width at half 
maximum (FWHM) of each line. 

To provide a better estimate of the FWHM of the H/3 
and [OIII] line profiles we also performed a fit of the spec- 
trum using five components model (Model 1): a power- 
law to account for the continuum emission, four Gaussian- 
shaped emission lines to account for the H/3 (both narrow 
and broad component) and the [OIII] 4959A, [OIII] 5007A 
lines. FWHM and velocity offset of the narrow lines are 
forced to be the same. Results are shown in Fig. 5 (upper 
panel). To account for residuals in the 4830 - 4900A range 
(i.e. the asymmetry in the broad H/3 line) we ran another 
fit by adding another emission line component named H/3* 
(Model 2, Fig. 5 lower panel). 

4.1 Results 

A qualitative comparison of the spectra by INT and SDSS 
taken ~9 years apart (Tab. 2) shows that the contirmum 
has changed both in intensity (by a factor ^5) and in slope. 
The intensity and shape of the Balmer line seem to be the 

^ http://www.sdss.org/dr7. Data from DR8 does not differ sig- 
nificantly from DR7. 



same in both spectra. In §3 we also analyze the broad- 
band SED (lower panel of Fig. 2) of the source, consider- 
ing also the photometric data from SDSS, taken ~1 year 
before the SDSS spectroscopic observations. Again, we see 
that a change in intensity of the optical emission of a fac- 
tor ~10 occurred in just 1 year. Also, a change in slope 
occurred, as indicated by the photometric data in the 5 
available bands. Thus, we conclude that B2 0954-I-25A has 
been observed in at least three emission states at optical 
wavelengths, which differs in intensity: from the dimmer one 
(SDSS/photometric), through an intermediate one (INT), to 
the brighter one (SDSS/spectroscopic). The spectral slope is 
correlated to the change of state. We will refer to such states 
as the disk dominated one, the intermediate state and the 
jet dominated state respectively, and discuss their properties 
in §3.3. 

A quantitative spectral analysis can only be performed 
on the SDSS spectrum. Fig. 4 clearly shows that our value 
of a=0. 70747 locates the peak of the [OIII] Une exactly at 
in velocity space, while the Mgll peak is rcdshiftcd by ~300 
km s~^ and the H/3 peak is blueshifted by ~100 km s~^. 
Also, the H/3 profile shows a pronounced asymmetry on the 
red wing. Rough estimates of the FWHMs are: 3050 km s~^ 
for the Mgll line, 2600 km s"^ for the H/3 line and ~480 km 
s"^ for the [OIII] line. 

To provide better estimates of luminosities, widths and 
offsets of emission lines, and eventually highlight some pecu- 
liarities in the spectrum, wc performed a fit with the models 
described above. The spectrum in the 4800 - 5100A region, 
as well as the fitting model and the single components, are 
shown in Fig. 5 (Model 1 in upper panel, Model 2 in lower 
panel). Results are shown in Tab. 5. 

Model 1: The resulting FWHM for the H/3 line is 
(2830 ± 220) km s"\ while the FWHM of the narrow lines 
is (430 ± 24) km s~^, grossly consistent with the rough esti- 
mates obtained with Fig. 4. The model fits reasonably well 
the spectrum, and the value is rather small. Also, the lu- 
minosity of the continuum at 5100A and the spectral slope, 
as well as the sum of luminosities from broad and narrow 
H/3 lines and [OIII] luminosities, are consistent with the fit 
performed in SIO (Tab. 5), thus our model does not suffer 
from lack of iron line modeling (which is present in their fit). 
The value of FWHM and velocity offset for both the broad 
and narrow components of H/3 aie however, quite different. 
The strong discrepancies in FWHM between our results and 
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Figure 3. Optical spectrum of B2 0954+25A observed by SDSS in 2006 (red line, flux scale on the left axis) and with the Isaac 
Newton Telescope (INT) in 1987 (black line, flux scale on the right axis). The range of values covered by both axes is exactly the 



same (~ 0.5 X 10 erg s 



^) so that spectra can be directly compared. SDSS spectrum have been de— reddened using a Milky 



Way— like extinction curve (Pei 1992) and Ab = 0.158 mag, transformed to the rest frame using our redshift estimate (z=0. 70747, see 
text) and rebinned by a factor of 2. The INT spectrum has been derived from the plot given in Jackson & Browne (1991, their Fig. 2) 
and transformed to rest frame. The "S" letter denote a telluric line. Some of the most important emission lines are highlighted. The 
spectra show that the continuum has changed both in intensity (by a factor ~5) and in slope. The intensity and shape of the Balmer 
lines seem to be the same in both spectra. 
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Figure 4. Proflles of Mgll (2800A), (4863A) and [OIII] (5007A) emission lines in velocity offset space. Spectrum has been continuum- 
subtracted, the continuum level has been estimated by eye and all proflles have been normalized to have maximum of 1. Finally, proflles 
have been smoothed by a ~4 A boxcar average to provide a clearer view. Dashed lines allow to roughly measure full width at half 
maximum (FWHM) of each line. 



© 2011 RAS, MNRAS 000, 1-14 



8 G. Calderone et al. 
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Figure 5. Spectrum of B2 0954+25A and fitted model in the range 4810—5100 A. All quantities are quoted in rest frame. Upper panel: 
spectrum (black points, la error bars) with fitted model (red points) and model components: broad lines (red solid line), narrow lines 
(blue solid lines) and power— law continuum (green solid line). The broad and narrow lines components are shifted upward for clarity. The 
five black dotted lines are the rest line centers of (from left to right) H7, H/3, "line 2" (see text), [OIII] 4959 and [OIII] 5007 respectively. 
Red and blue dotted lines are the rest line centers of all broad and narrow considered lines. Lower panel: residual in units of one standard 
deviation. Red solid line is the cumulative reduced chi-squared (values on right axis). 
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those from SIO are probably due to the fact that wc had the; 
opportunity to carefully check the results of the fit, while in 
SIO an automatic fitting pipeline had to be used to handle 
the many spectra contained in the catalog. 

Model 2: The residuals of the fit in Model 1 (Fig. 5, 
upper panel) look random everywhere, except in two re- 
gions: 4840-4900A and 4930-4945A, in which the residu- 
als are quite "coherent" , suggesting the presence of further 
components. Although the model is far from being rejected 
{Xied < 1) these residuals are due to the asymmetry in 
the broad H/3 profile discussed above, and to an additional 
emission line. Thus, we ran again the fit with a new broad 
H/3-like component, named H/3* in Tab. 5 (actually we are 
not interested in the additional line at 4940A). Results are 
shown in Fig. 5 (lower panel) and Tab. 5. The extra com- 
ponent accounts well for the asymmetry of the H/3 profile, 
and the residuals in the region now look random. The same 
consideration as for Model 1 applies for the luminosity and 
slope of the continuum, and the luminosity of the H^ and 
[OHI] lines. Clearly, the putative broad H^ line is now much 
narrower (~1500 km s~^). Furthermore, the new H^S* com- 
ponent has the same FWHM and velocity offset as the broad 
H/3 values given in SIO. This provides a possible explanation 
for the discrepancies between our results and those from SIO. 

In summary, our Model 2 seems to be the best descrip- 
tion of the spectrum of B2 0954-I-25A in the 4810 - 5100 A 
range. 



5 THE LINE PROFILE 

The FWHM value of the H/3 line profile reported in litera- 
ture is 65 A (Jackson & Browne 1991). As shown in Tab. 

6 the FWHM found with both Model 1 and 2 of our fit is 
significantly lower (46 A and 24 A respectively). Since the 
H/3 profiles observed with INT and SDSS (Fig. 3) overlap 
perfectly, we can exclude a variation in the line profile, and 
conclude that previous value was likely overestimated. More- 
over, the value reported in SIO is perfectly compatible with 
our estimate in Model 1 (48 A, Tab. 6). 

Fig. 5 show that the broad H/3 emission line can hardly 
be modeled with a single Gaussian profile, nor with any 
symmetric line profile such as the often quoted logarith- 
mic profile (Blumenthal & Mathews 1975), since the red 
wing is highly asymmetric. Iron emission lines in the range 
4840-4900 A are usually much weaker than the Balmer lines 
(Capriotti et al. 1979), thus such profile is hardly the result 
of a blending of diS'crcnt lines. Iron lines may be present 
but their intensity is negligible, indeed our fitting procedure 
provide the same results for line luminosities and width of 
narrow lines as reported in SIO (§4.1). Asymmetric H/3 line 
profiles are not new in Type 1 AGN spectra: both red and 
blue asymmetries has been observed in all kind of AGNs (e.g. 
Capriotti et al. 1979, 1980; Peterson et al. 1987; Stirpe 1990; 
Romano et al. 1996; Veron-Cetty et al. 2001), as well as pro- 
file variability (Osterbrock & Phillips 1977; Peterson 1987; 
Stirpe et al. 1988). Several kinematic models have been pro- 
posed to explain such features in the broad Balmer spectra of 
AGN, including radial inflow or outflow, different geometries 
and multiple components BLR (e.g. Peterson 1987; Popovic 
et al. 2004; Romano et al. 1996; Zhu et al. 2009), but no 
one has yet reached a general consensus. This hints for the 



presence of non virialized components in the BLR region. 
Thus, measures of FWHM performed on such profiles may 
not be directly related to the underlying black hole mass. 
This is the reason why wc used the FWHM of just the main 
H/3 line to estimate the virial black hole mass, neglecting the 
additional H^* component which account for the asymme- 
try. Doing so, the FWHM of H^ turns out to be about half 
the FWHM of Mgll (Fig. 4). Such differences are common 
in sources with small values of FWHM(H/3): in the SIO cat- 
alog, the Mgll line is wider than H/3 for ~84% of the sources 
with FWHM(H/3)< 2000, and only 22% of the sources with 
FWHM(H/3)> 4000. 

Of course the reasoning can be inverted: the "real" H/3 
line profile may be much wider, but some intervening gas 
could absorb radiation only on the blue side. Although we 
have no evidence to exclude this hypothesis, the range of 
values for the black hole mass found in §5.1 marginally sup- 
ports the hypothesis that the real virialized line profile is the 
one modeled with the main H^S component. This conclusion 
may have important consequences on the classification of 
B2 0954-I-25A. Should the additional H/3* component vanish 
in future observations, the FWHM of the remaining, virial- 
ized, H/3 component would be ~1500 km s^^, and the source 
would become a powerful 7-RL-NLSI. (Foschini et al. 2009; 
Abdo et al. 2009c; Calderone et al. 2011). Such variations 
in line profiles have already been observed, on timescales of 
tens of years in NGC 5548 (Peterson 1987). 

5.1 Spectroscopically derived quantities 

Spectroscopic measurements, such as those derived in §4.1, 
are used to infer some physical jiropcrties of the source, 
such as the total accretion luminosity Lboi and the virial 
black hole mass. The former is obtained by scaling the lu- 
minosity of one emission line (e.g. H,/3) to obtain the total 
BLR luminosity, Lblr, through a template quasar spectrum 
(Francis et al. 1991). Lboi can then be computed by divid- 
ing Lblr by the covering factor, assumed to be 0.1. Black 
hole masses are usually estimated assuming virialized mo- 
tion of BLR clouds. The BLR radius has been measured by 
means of the reverberation mapping technique (RM, Pe- 
terson 1993) for a few tens of objects, and a correlation be- 
tween the continuum optical luminosity and the BLR size, 
has been calibrated [R-L relation, Kaspi et al. 2005; Vester- 
gaard & Peterson 2006; Bentz et al. 2009), so that wc can 
infer Rb^r size without the (very time-consuming) RM ob- 
servational campaigns. Thus, it is in principle possible to 
estimate black hole masses using single-epoch optical spec- 
troscopy (such as the SDSS spectrum of B2 0954+25A used 
in §4) and a calibrated virial mass scaling relationship (such 
as those reported in Vestergaard & Peterson (2006)). There 
are, however, several caveats, such as the geometry and ori- 
entation of the BLR (Decarii et al. 2011), the role of radi- 
ation pressure (Marconi et al. 2008, 2009) and the way line 
profiles should be analyzed to infer the orbital velocity of 
the clouds. Indeed, the H/J broad line is often asymmetric, 
thus the broadening cannot be simply interpreted as being 
due to orbital bulk motion. In these cases is not clear what 
measure should be taken for the line width. Note also that 
RM studies have been carried out mainly on radio-quiet 
sources, i.e. sources in which synchrotron emission from a jet 
at optical/IR wavelengths is negligible with respect to the 
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Table 5. Quantities derived from fit parameters and corresponding quantities in the Shen et al. (2010) catalog. 



Parameter 


Model 1 


Model 2 


SIO 


Units 




Lum. (H/3broad) 


26.3 ± 2.2 


12.9 ± 6.3 


23 ± 13 


[10''2 erg s- 




Lum. (H/3*) 


— 


12.1 ± 5.8 


— 


[10*2 erg s- 




Lum. (H/3narrow) 


0.93 ± 0.58 


0.50 ± 0.92 


8 ± 42 


[10''2 erg s- 




Lum. ([OIII], 5007) 


11.07 ± 0.59 


11.01 ± 0.58 


12.2 ± 1.9 


[10''2 erg s- 




FWHM (H/3b,oad) 


2830 ± 220 


1470 ± 320 


1870 ± 600 


[km s~^] 




FWHM (11;'^*) 




1790 ± 600 




[km s-1] 




FWllAi (llJuarrow) 


l:« ± 2 i 


i2« ± 2 i 


1200 ± 100 


[km.s i; 




VofI (H/3broad) 


344 ± 75 


-(350 ± 270) 


-(1250 ± 630) 


[km s^^] 




VotF (H/3*) 




1220 ± 400 




[km s"-*-] 




Voff (H/3narrow) 


-(0 ± 10) 


1 ± 10 


37 ± 200 


[km s"-*-] 




ALa (5100A) 


220.88 ± 0.80 


220.59 ± 0.78 


220.73 ± 0.74 


[10** erg s" 




Continuum index 


-(0.75 ± 0.13) 


-(0.85 ± 0.11) 


-(0.709 ± 0.018) 







emission from the so called Big Blue Bump (BBB). Thus, 
the R—L relations are reliable only when the BBB is clearly 

visible (Wu ct al. 2004). Contribution from host starlight 
is also supposed to be an issue (e.g. Bentz ct al. 2009), al- 
though only in low luminosity AGN (log(ALA,5100A)<44.5, 
Shen et al. 2010), therefore not in the case of B2 0954-I-25A. 
Finally, the \a statistical dispersion of virial mass scaling 
relationships against RM based masses is estimated to be 
~0.5 dex. The RM masses are in turn scattered by the same 
factor around the M-a relation (Onken et al. 2004). The 
absolute accuracy is therefore estimated to be of the order 
of 0.7 dex (Vcstcrgaard & Peterson 2006). 

We will consider virial mass estimates computed us- 
ing both Model 1 and 2. Also, the SDSS spectrum of B2 
0954+25A has been taken when the source was in the jet 
dominated state, therefore we cannot use the continuum lu- 
minosity estimated from spectral analysis to infer the BLR 
size. Instead, we will use optical data of the disk dominated 
state, as observed photometrically by SDSS. 

5.1.1 Results 

To compute Lblr, and the disk bolometric luminosity Lboi, 
we use the composite quasar spectrum given in Francis et al. 
(1991), including the Ha contribution as in Celotti et al. 
(1997). The broad H/3 line contributes to the entire Lblr 
with a weight 22/555 (the Lya contributes for 100/555). 
Using the luminosity of the H/3 complex found with Model 
2 fitting (broad H/3 line + H/3* component, log(L(H/3)/crg 
s~^) = 43.4), and setting Lboi = IOLblr, we estimate 
log(Lboi/erg s~^) = 45.8. This value agrees with the one 
derived in §3.2 through SED modeling: log(Lboi/erg s~^) 
= 46. Instead, SIO reports log(Lboi/erg s^^) = 47, because 
SIO used an observation made when the source was in the 
jet dominated state. 

Tab. 6 reports the black hole mass estimates for B2 
0954+25A obtained using different methods and calibra- 
tions. For all estimates, except the last two, we used the 
virial mass scaling relations. 

The first two estimates are computed using the FWHM 
estimates (of the broad component and the whole profile 
respectively) of the Mgll emission line provided by SIO, the 
continuum luminosity derived from SDSS photometry (when 
the source was in the disk dominated state) of log{i'L„/eTg 
s~^) = 45.7 and the same calibrations used in SIO: Shen 
et al. (2010) and Vestergaard & Osmer (2009) respectively. 



The mass estimates are very similar, so we will refer to them 
as the Mgll virial mass estimates (log M/Mq ^ 9). 

The next two estimates are computed using the FWHM 
of the H/3 emission line fitted with a single Gaussian. Values 
for FWHM are provided by SIO and by our Model 1 respec- 
tively, the continuum luminosity is derived from SDSS pho- 
tometry (disk dominated state): log(:/L^/erg s~^) = 45.3, 
and the calibration are from Vestergaard & Peterson (2006, 
same as SIO) and Bentz et al. (2009), respectively. This 
pair of mass estimates are also very similar, so we will 
refer to them as the single H^ virial mass estimates (log 
M/Mq ~ 8.5). 

Next estimate is computed using the FWHM of the H/3 
emission line from Model 2, neglecting the additional H/3* 
component, and the calibration by Bentz et al. (2009). 

Finally, the last two rows of Tab. 6 report mass esti- 
mates obtained without using the virial method. For the 
first one we used the FWHM of the [OIII] line (provided by 
Model 2 fitting) as a proxy for the stellar velocity dispersion, 
and the calibration of the M — a* relation given in Tremaine 
et al. (2002, their Eq. 20). In the last one we simply report 
the mass estimate obtained in §3.2 through SED modeling. 
Notice that the last three methods yield very similar results 
(log M/Mq ~ 8), even if they have been derived using dif- 
ferent and independent methods. 

5.1.2 Discussion 

Our black hole mass estimates are considerably lower than 
values found in literature (§2). This is likely a consequence 
of the different values of the H,0 line profile's width, as dis- 
cussed in 5. Also, the mass given in SIO is probably overes- 
timated since the continuum luminosity has been measured 
while the source was in the jet dominated state. Therefore, 
we will consider only our mass estimates given in Tab. 6. 

In particular, the estimate obtained with Model 2 (ne- 
glecting the H/3* component), and the Af — ct* relation 
(log(M/MQ) ~8) are instead quite independent, thus they 
provide greater confidence about the goodness of the mass 
estimate. Modeling the disk dominated state with a standard 
(i.e. Shakura & Sunyaev 1973) disk (§3.3) yields estimates 
of the black hole mass between (1.5 - 5) x IO^Mq. Smaller 
masses are excluded because the disk would become super- 
Eddington, larger masses are excluded because, beside pro- 
viding a worse fit to the near-IR and optical data of the 
"disk-dominated" state, they require a smaller ionizing lu- 
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Table 6. Black hole mass estimates using different methods and calibrations. Continuum luminosity at 
3OOOA and 5IOOA (rest frame) used in virial estimates are log{!^L;^/erg s~^) = 45.7 and 45.3 for Mgll and 
H/3 respectively (derived from SDSS photometry). Col. [1]: method used; Col. [2]: reference for the FWHM 
of the emission line; Col. [3]: FWHM of the broad emission line; Col. [4]: black hole mass estimate; Col. [5]: 
Eddington ratio, assuming a bolometric luminosity log(Li,oi/erg s~^) = 45.95 (see §3.2 and Tab. 3); Col. [6]: 
calibration references: (a) Shen et al. (2010), (b) Vestergaard & Osmer (2009), (c) Vestergaard & Peterson 
(2006), (d) Bentz et al. (2009), (e) Tremaine et al. (2002). 



Method 


FWHM ref. 


FWHM [km s-\ (A)] 


log M/Mq 


Edd. ratio 


Calib. ref. 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


Virial, Mgll 


SIO, broad comp. 


3980 ( 65 ) 


9.0 


0.1 


(a) 


SIO, whole profile 


3390 ( 55 ) 


8.8 


0.1 


(b) 


Virial, H/3, single comp. 


SIO 

Model 1 


2970 ( 48 ) 
2830 ( 46 ) 


8.6 
8.5 


0.2 
0.3 


(c) 
(d) 


Virial, H^, H^* 


Model 2 


1470 ( 24 ) 


7.9 


1.1 


(d) 


M-0-, ([OIII]) 


Model 2 


430 ( 7 ) 


8.0 


0.9 


(e) 


SED modeling 






8.2 


0.6 
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Figure 6. Comparison of the SED of B2 0954+25A with the SED 
of 3C 273 (top panel) and PMN J0948+0022 (bottom). 

minosity, that becomes incompatible with the luminosity of 
the broad line region. We therefore conclude that the black 
hole mass of B2 0954+25 A is in the range log(M/M0) = 
8-8.5. 



0954+25A, although slightly steeper. At the lowest observed 
frequencies the two SED are markedly different: whereas the 
74-365 MHz radio spectral index of 3C 273 become steep 
(q ~ —0.5), the spectral index of B2 0954+25A remains 
flat (a ~ —0.2). As discussed in §3.1, this behaviour can 
be explained with a different amount of relativistic beaming 
which boosts the jet emission in B2 0954+25A. At shorter 
wavelengths 3C 273 is ~1 order of magnitude more power- 
ful than B2 0954+25A. If emission in both sources occur at 
approximately the same Eddington ratio (0.4-0.6, Tab. 3) 
as suggested by SED modeling, the luminosity ratio at opti- 
cal/UV wavelengths (~ 3) provides an estimate for the mass 
ratio. Since log(M/M0) of 3C 273 is 8.9, for B2 0954+25A 
we expect a mass of ~8.4. This is yet another indication that 
the black hole mass of B2 0954+25A should be significantly 
smaller than log(M/MQ) = 9. 

Also, the SED of PMN J0948+0022 is very similar to 
that of B2 0954+25A (in the disk dominated state), the only 
difference being at 7-rays (where PMN J0948+0022 showed 
an exceptional flaring episode (Foschini et al. 2011)). Pa- 
rameters from the SED modeling and jet powers are quite 
similar (Tab. 3 and 4). This similarities support our spec- 
ulation on the possible classification of B2 0954-I-25A as a 
7-NLSI. 

Ifowever, this similarity cannot be pushed too far, since 
B2 0954-(-25A lacks the strong iron emission that usu- 
ally characterizes NLSl. Also, X-ray photon index of B2 
0954-f25Ais F = 1.74 (§3), while that of NLSl is typically 
F > 2.5 (Boiler et al. 1996). However, consider that also 
the other NLSl detected by Fermi have a flat X-ray spec- 
trum (Abdo et al. 2009b, c) likely because the flat X-ray jet 
emission dominates. 



6 COMPARISON WITH 3C 273 AND PMN 
J0948+0022 

In Fig. 6 we compare the broad band SED of B2 0954-I-25A 
to that of 3C 273 (upper panel) and PMN J0948-I-0022 
(lower panel) . Both sources have been analyzed in Ghisellini 
et al. (2010) and model parameters are reported in Tab. 3 
and 4. 

The radio spectrum of the prototypical blazar 3C 273 
shows a flat spectrum comparable in luminosity with B2 



7 TYPICAL FERMI BLAZAR OR 7-NLSI? 

As discussed in §3.2 and §6, B2 0954+25A is a typical Fermi 
blazar. However, the likely small black hole mass of B2 
0954-I-25A and the SED resemblance with prototypical 7- 
NLSl PMN J0948-I-0022 suggest some similarities between 
the two sources. B2 0954+25A does not meet the commonly 
adopted classification criterion for NLSl, namely a FWHM 
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of broad H/3 <2000 km s'^ (Osterbrock & Pogge 1985).* 
However, the physical nature of such empirical threshold 
has been often questioned (e.g. Goodrich 1989; Veron-Cetty 
et al. 2001) since all observational properties show a continu- 
ous transition at FWHM(H/3) ^2000 km s~^, i.e. properties 
of NLSl and BLSl sources are smoothly joined. 

Alternative criteria to distinguish among BLSl and 
NLSl sources have been put forward. Sulentic et al. (2000) 
proposed to increase the dividing threshold to 4000 km s^^, 
to select radio-quiet sources with low mass and high accre- 
tion rate. Veron-Cetty et al. (2001) suggested to consider the 
strength of Fell emission relative to H/3 as a possible tracer 
of the Eddington ratio. Netzor & Trakhtenbrot (2007) pro- 
posed to classify as narrow-line AGN (NLAGN) all sources 
exceeding an Eddington ratio of 0.25 (regardless of black 
hole mass). 

A small black hole mass is likely a characterizing prop- 
erty of NLSl sources (Grupe & Mathur 2004). We would like 
to stress that if asymmetry in broad Balmer line profiles are 
due to non-virialized components in the BLR, then only the 
virialized one should be considered to estimate the black 
hole mass, and thus the NLSl classification. In this case B2 
0954+25A would be classified as a powerful 7-NLSI, just 
like PMN J0948-h0022. 

Since the overall properties of B2 0954-I-25A are similar 
to both classes of Fermi blazars and 7-NLSI, we conclude 
that it is a one of those objects which smoothly joins the 
population of narrow-line and broad-line AGN. 



8 SUMMARY AND CONCLUSIONS 

In this paper we carried out an extensive data analysis on 
the Fermi blazar B2 0954+25A using archival data, to point 
out some peculiarity and try to provide coherent picture of 
its physical properties. The main conclusions are as follows: 

• (§2, §3) the source B2 0954-I-25A is a strong 7-ray emit- 
ter, detected in both the lyr and 2yr Fermi/LAT point 

source catalogs (Abdo et al. 2010, 2011). We found no sig- 
nificant variability of the flux nor of the spectral index in 
the 0.1-100 GeV band in the first 34 months of data from 
Fermi/LAT; 

• (§3.1) its radio spectrum is "flat" down to very low 
frequencies (74 MHz, Cohen et al. 2007), i.e. the optically 
thick synchrotron emission from the jet is dominating over 
the optically thin, isotropic emission from extended region. 
This observation indicates that the viewing angle must be 
small (3-6 degrees). The corresponding de-projected size of 
the jet is in the range 0.5-1 Mpc, the size of a giant radio 
lobe. The isotropic and the beamed components should have 
equal luminosities at ~10 MHz; 

• (§3.1, Fig. 1) the "flatness" of the radio spectrum at low 
frequencies observed in B2 0954+25A is typical of Fermi 
detected blazars (2LAC sample), i.e. well aligned blazars. 
Radio selected blazars show steeper spectra on average (e.g. 
the CRATES sample). 



* The further criterion on flux ratio [OIIIJ/H/9 <3 is almost al- 
ways automatically satisfied when a broad H/3 component is de- 



• (§4.1, Fig. 3 and 2) B2 0954-I-25A has been observed in 
at least three different optical states: from a disk-dominated 
state, to a jet dominated one in which the synchrotron emis- 
sion overwhelms the emission from the disk, via an interme- 
diate state. The two extreme states differ by nearly one order 
of magnitude in luminosity. Despite the change in contin- 
uum luminosity, the broad hydrogen Balmer lines maintain 
the same luminosity and line proflle. 

• (§4.1, Fig. 5, Tab. 5) the previously known estimate for 
the FWHM of the broad H/3 line proflle (65 A rest frame, 
corresponding to ~4000 km s~^) from Jackson & Browne 
(1991) is probably overestimated. Also, the value of 1870 
km s^^ reported in Shcn et al. (2010) is probably due to a 
poor model fit. A more likely value is 46 A, corresponding 
to ~2800 km s~^. The broad H/J line shows a pronounced 
asymmetry on the red wing, which require a second Gaussian 
component (H/3*) of similar FWHM, to be modeled. The 
velocity offset of the H/3* is ~1200 km s"^ 

• (§4.1, Fig. 4) [OIII] 5007A and Mgll lines are quite 
symmetric. The FWHM of the Mgll line is ~ 3050 km s"\ 
and the peak is redshifted by ~300 km s~^ with respect to 
the redshift identified by the [OIH] line. 

• (§5.1, §5.1.2) the most likely value of the black hole mass 
is M ~ (1 — 3) X IO^Mq, as suggested by several independent 
methods. In particular, for the mass scaling relations we 
considered only the symmetric component of the H/3 profile, 
and the continuum luminosity in the disk dominated state 

• (§3.2, Fig. 2, Tab. 3, §6) modeling the broad band spec- 
tral energy distribution, we derive physical parameters quite 
typical of other Fermi blazars. 

• (§6) We suggest to classify B2 0954-I-25A as a transi- 
tion object between the class of FSRQ and 7-NLSI, since 
it shows characteristic features of both classes (namely, the 
blazar appearance and the similarity with the SED of PMN 
J0948-F0022). 



APPENDIX A: APPENDIX: SOME DETAILS 
ON THE MODELING 

We use the model described in detail in Ghisellini & Tavec- 
chio (2009). The emitting region is assumed spherical, at a 
distance Rdiss from the black hole, of size R = ipRnsB (with 
ip — 0.1) from the black hole of mass M, moving with a bulk 
Lorcntz factor F. The bolometric luminosity of the accretion 
disk is id- 

The energy particle distribution N{'y) [cm~^] is cal- 
culated solving the continuity equation where particle in- 
jection, radiative cooling and pair production (via the 7- 
7 — i> process), arc taken into account. The created pairs 
contribute to the emission. The injection function Q{'y) 
[cm~^ s~^] is assumed to be a smoothly joining broken 
power-law, with a slope Q{'y) oc 7"''-'^ and 7~^^ below and 
above a break energy 7b: 



Q(7) = Qo 



l + (7/7b)-''i+"^ 



(Al) 



tected (e.g. Zhou et al. 2006). 



The total power injected into the source in the form of rel- 
ativistic electrons is P( = rricC^V J Q{'l)^d'y, where V = 
{4Tr/3)R^ is the volume of the emitting region. 

The injection process lasts for a light crossing time R/c, 
and we calculate N{'y) at this time. This assumption comes 
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from the fact that even if injection lasted longer, adiabatic 
losses caused by the expansion of the source (which is trav- 
eling while emitting) and the corresponding decrease of the 
magnetic field would make the observed flux to decrease. 
Therefore the calculated spectra correspond to the maxi- 
mum of a flaring episode. 

Above and below the accretion disk, in its inner parts, 
there is an X-ray emitting corona of luminosity Lx (it is 
fixed at a level of 30% of Ld). Its spectrum is a power law 
of energy index ax = 1 ending with a exponential cut at 
Ec =150 keV. The specific energy density (i.e. as a function 
of frequency) of the disk and the corona are calculated in the 
comoving frame of the emitting blob, and used to properly 
calculate the resulting External inverse Compton spectrum. 
The BLR is assumed to be a thin spherical shell, of radius 
JJblr = W^^L\^^^ cm. We consider also the presence of a 
IR torus, at larger distances. The internally produced syn- 
chrotron emission is used to calculate the synchrotron self 
Compton (SSC) flux. Table 3 lists the adopted parameters. 

Table 4 lists the power carried by the jet in the form of 
radiation (Pr), magnetic field (Pb), emitting electrons (Pe, 
no cold electron component is assumed) and cold protons 
(Pp, assuming one proton per emitting electron). All the 
powers are calculated as 



(A2) 



where U'i is the energy density of the i component, as mea- 
sured in the comoving frame. 

The power carried in the form of the produced ra- 
diation, Pr — ttR^F^ PcU^^, can be re-written as [using 
[7;,d = i7(4rfc)]: 



Pr = L 



4J4 



4<52 



(A3) 



where L is the total observed non-thermal luminosity (L' 
is in the comoving frame) and f/^ad is the radiation energy 
density produced by the jet (i.e. excluding the external com- 
ponents). The last equality assumes 0^ ~ l/F. 

When calculating Pe (the jet power in bulk motion of 
emitting electrons) we include their average energy, i.e. Ui = 
ne(7)meC^. 
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